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"j^ , Abstract. We study the phase diagram of the 2D Ji-J[-J2 spin-1/2 Heisenberg 

model by means of the coupled cluster method. The effect of the coupling J[ on the 
Neel and stripe states is investigated. We find that the quantum critical points for 
the Neel and stripe phases increase as the coupling strength J[ is increased, and an 
Q ' intermediate phase emerges above the region at J[ « 0.6 when Ji — 1. We find 

O i indications for a quantum triple point at J[ w 0.60 ± 0.03, J2 ~ 0.33 ± 0.02 for Ji = 1. 
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1. Introduction 



O 

(N 

^ ! The spin-1/2 Heisenberg antiferromagnet has been much studied in recent years due 
to the discovery or successful syntheses of such new magnetic materials as the layered- 



! oxide high-temperature superconductors. Also of much interest has been the interplay 
^ ' between frustration and quantum fluctuations in 2D quantum spin systems that can 
lead to quantum phase transitions between magnetically ordered semiclassical and novel 
^ ! quantum paramagnetic ground-state phases, see, e.g., Refs. [1,2]. 
5^ i For example, the frustrated 2D antiferromagnetic J1-J2 niodel with nearest- 

neighbour (Ji) and next-nearest-neighbour (J2) bonds has attracted much attention 
both theoretically (see, e.g., Refs. [3-10] and references cited therein) and 
experimentally [11, 12]. It is now well accepted that the model exhibits two phases 
displaying semiclassical magnetic long-range order (LRO) at small and at large J2, 
separated by an intermediate quantum paramagnetic phase without magnetic LRO in 
the parameter region Jci < J2 < Jc2 where Jc^ ~ 0.4Ji and ~ O.6J1. The ground 
state (GS) for J2 < Jci exhibits Neel magnetic LRO, whereas for J2 > it exhibits 
collinear stripe LRO. 

In real systems deviations from the ideal 2D J1-J2 model such as spin 
anisotropies [8,13] or interlayer coupling may be relevant [10]. In addition, a 3D version 
of the J1-J2 model [14] has also been considered. 
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An interesting generalization of the pure J1-J2 model has been introduced recently 
by Nersesyan and Tsvelik [15]. They consider a spatially anisotropic spin-1/2 2D Ji- 
J[-J2 model, where the nearest-neighbour bonds have different strengths Ji and J[ in 
the X and y directions. This model has been further studied by other groups using 
the exact diagonalization (ED) of small lattice samples of < 36 sites [16], and the 
continuum limit of the model [17]. Both groups support the prediction of the resonating 
valence bond state by Nersesyan and Tsvelik [15] for J2 = 0.5J( ^ Ji, and the limit of 
small interchain coupling extends along a curve nearly coincident with the line where 
the energy is maximum. The model has also been studied by Moukouri [18] using a 
two-step density-matrix renormalization group approach. 

Our aim here is to further the study of this model by using the coupled cluster 
method (CCM). The CCM (see, e.g., Refs. [19-21] and references cited therein) is one of 
the most powerful and universally applicable techniques of quantum many-body theory. 
It has been applied successfully to calculate with high accuracy the ground- and excited- 
state properties of many lattice quantum spin systems (see, e.g., Refs. [10,21-26] and 
references cited therein). A particularly important result from our calculations is the 
indicated existence of a quantum triple point (QTP) at nonzero (positive) values of Ji, 
J[ and J2. 



2. The model 



The Ji-J(-J2 model is a spin-1/2 Heisenberg model on a square lattice with three 
kinds of exchange bonds, with strength Ji along the row direction, J[ along the column 
direction, and J2 along the diagonals, as shown in figure [H All exchanges are assumed 
positive here, and we set Ji = 1. The Hamiltonian of the model is described by 

H = Ji ^ Sj^; ■ Sj+i^/ -I- J[ ^ Sj^/ ■ Sj^;+i 
i,l i,l 

+ J2 "^{^l ■ + Si+i,i ■ Si,i+i), (1) 

i,l 

where the index (i, /) labels the x (row) and y (column) components of the lattice sites. 

This model has two types of classical GS, namely, the Neel (vr, vr) state and stripe 
states (columnar stripe (vr, 0) and row stripe (0,7r)), the spin orientations of which are 
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Figure 1. (a) Ji-J[-J2 model; — Ji; - - J(; • — ■ J2; (b) Neel state, (c) stripe state - 
columnar and (d) stripe state - row. Arrows in (b), (c) and (d) represent spins situated 
on the sites of the square lattice (indicated by • in (a)). 
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shown in figures[T](b,c,d) respectively. There is clearly a symmetry under the interchange 
of rows and columns, Ji ^ J[, which implies that we need only consider the range of 
parameters with J[ < Ji. The ground-state (gs) energies of the various classical states 
are given by 



N 4 
1 



columnar 



4 

E^' 1 



row 



Jl- J(+2J2), 
(-Ji + J(-2J2), 

^ ^(J,-J[-2J2). (2) 

We take Ji = 1 and J[ < 1. Clearly, from ([2]), the classical GS is then either the 
Neel state (if J[ > 2J2) or the columnar stripe state (if J[ < 2J2). Hence, the (first- 
order) classical phase transition between the Neel and columnar stripe states occurs at 



3. The coupled cluster method 

The CCM formalism is now briefly described (and see Refs. [19-26] for further details). 
The starting point for the CCM calculation is to select a normalized model state |$). 
It is often convenient to take the classical ground state as the model state for spin 
systems. Hence our model states are the Neel state and the columnar stripe state. In 
order to treat each site identically, we perform a mathematical rotation of the local 
axes of the spins such that all spins in the reference state align along the negative z- 
axis. The Schrodinger ground-state ket and bra CCM equations are -f^l^E') = -E|\E') and 
(\E'|/f = £'(\E'| respectively. The CCM employs the exponential ansatz, =e'^|<l>) and 
(\E'| = (^liSe""^. The correlation operator S is expressed as S* = Xl/^o'^^^/ 
counterpart is 5 = ^ + J2i^o^i^7 ■ The multispin creation operators Cf = {CyY, with 
Cq = 1, are written as = sj^s^ ■ ■ ■ s^. The gs energy is E = ($|e~'^ife'^|$); and the 
staggered magnetization M in the rotated spin coordinates is M = — -^(^| ^jl"^)- 
The ket- and bra-state correlation coefficients {Si,Sj) are calculated by requiring 
the expectation value H = (\^|iJ|\E') to be a minimum with respect to all parameters 
{Si, Si) such that ($|C7e-^i/e^|$) = and {^\S{e-'^He^ - ^o)C/|$) = ; V/ ^ 0. 



4. Approximation scheme 

The CCM formalism is exact if all spin configurations are included in the S and S 
operators. In practice, however, truncations are needed. As in much of our previous 
work we employ here the localized LSUBn scheme, [6,10,21-26] in which all possible 
multi-spin-flip correlations over different locales on the lattice defined by n or fewer 
contiguous lattice sites are retained. The numbers of such fundamental configurations 
(viz., those that are distinct under the symmetries of the Hamiltonian and of the model 
state 1$)) that are retained for the Neel and stripe states of the current model in various 
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Table 1. Number of fundamental LSUBn configurations (fl f.c.) for the Neel and 
stripe states of the spin-1/2 Ji-J[-J2 model. 



Method 


tl 


f.c. 




Neel 


stripe 


LSUB2 


2 


1 


LSUB4 


13 


9 


LSUB6 


146 


106 


LSUB8 


2555 


1922 


LSUBIO 


59124 


45825 



LSUBra approximations are shown in Table [T] In order to solve the corresponding 
coupled sets of CCM bra- and ket- equations we use parallel computing [27]. 

5. Extrapolation scheme 

In practice one needs to extrapolate the raw LSUBn data to the n ^ oo limit. Based 
on previous experience [10, 23, 26] we use the following empirical scaling laws for the 
extrapolations of the gs energy per spin, 

E /N = flo + CLin^'^ + a2n~^ , (3) 

and the staggered magnetization of frustrated models, 

M = 6o + hn-'' , (4) 

where the exponent i/ is a fitting parameter. In order to fit to any fitting formula that 
contains n unknown parameters, one desirable rule is to have at least (n + l) data points 
(n + 1 rule) to obtain a robust and stable fit. In our results below the LSUBn results 
for n = {4, 6, 8, 10} are extrapolated. 

6. Results 

In figure [2] we show the gs energy per spin as a function of J2- For each value of J( 
two curves are shown, one (for smaller values of J2) using the Neel state, and the other 
(for larger values of J2) using the stripe state as CCM model state. Both sets of curves 
have the natural termination points [6,21,22] shown. For J(< 0.6 the two curves for a 
given value of J[ cross (or, in the limit, meet) very smoothly near their maxima, at a 
value of J2 slightly larger than the classical transition point of 0.5J(. This behaviour is 
indicative of a second-order quantum phase transition, by contrast with the first-order 
classical transition from ([2]). For J[ > 0.6 the curves no longer cross at a physical value 
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Figure 2. Extrapolated CCM LSUBn results for the gs energy per spin, E/N, for 
J[ = 0.2,0.4,0.6,0.8,1.0, using the Neel and stripe states of the s = 1/2 Ji-J[- 
J2 model. The LSUBn results are extrapolated in the limit n —^ 00 using the set 
n = {4, 6, 8, 10}. The NN exchange coupHng Ji = 1. 



(viz., -where the calculated staggered magnetization is positive), indicating the opening 
up of an intermediate quantum phase between the semiclassical Neel and stripe phases. 

In figure [3]^a) we show the equivalent staggered magnetization, M, for the Neel state, 
extrapolated using (jlj). We observe that the extrapolation scheme produces smooth and 
physically reasonable results, except for a very narrow anomalous "shoulder" region near 
the points where M vanishes for 0.6 ^J[^ 0.75. This critical regime is undoubtedly 
difficult to fit with the simple two-term scheme of (jl]). Our method for curing this 
problem and stabilizing the curves is to make efficient use of the information we obtain 
in (jlj) to extract the exponent z/, and then to use that value to infer the next term in the 
series. We find, very gratifyingly, that the value for u fitted to (jl]) turns out to be very 
close to 0.5 for all values of J[ and J2 except very close to the critical point. Therefore, 
we use the form of (I5l), 

M = Co + (ci + C2n-i) . (5) 

The use of ([5]) now removes the anomalous shoulder, as shown in figure [3]^b). Henceforth, 
in all of the results we discuss, we use (jS]) for the staggered magnetization. 

The raw LSUBn data for ?7.={2,4,6,8,10} and n={6,8,10} are also extrapolated. The 
results for n={2,4,6,8,10} and ?7,={4,6,8,10} are similar as they both obey the n+1 rule 
as mentioned above. This adds credence to the validity and stability of our results. The 
extrapolated curve for n = {6, 8, 10} has a very minor "shoulder" which is undoubtedly 
due to using only three data points to fit the three unknown terms, thus violating the 
n+1 rule. 

In figure H] we show our equivalent results for the staggered magnetization to those 
in figure [2] for the gs energy. We note the surprising result that M vanishes for both 
the quantum Neel and stripe phases at almost exactly the same critical value of J2, for 
a given J[, so long as J[ < 0.6. Conversely, for J[ > 0.6 there exists an intermediate 
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Figure 3. Extrapolated CCM LSUBn results for the gs staggered magnetization, M, 
for J( = 0.6,0.65,0.7 for the Neel state of the .s = 1/2 Ji-J[-J2 model, (a) Results 
using (gl), M = 6o + bin^" . (b) Resuhs using M = cq + n-°-^ (ci + C2n~^). The 
LSUBn results are extrapolated in the limit n ^ oo using the set n = {4,6,8,10}. 
The NN exchange coupling Ji — 1. 



region between the critical points at which M — for these two phases. The order 
parameters of both the Neel and the stripe phases vanish continuously both below and 
above the correspondingly indicated quantum triple point (QTP), as is again typical 
of second-order transitions. We note, however, that there exists some evidence from 
such other sources as ED calculations for the J1-J2 model (i.e., the present model with 
J[ = Ji) in 2D [5], that the transition between the stripe and intermediate phases is 
first-order, while others [28] have argued it may be close to second-order. We discuss 
these points further in Sec. [7] below. 

We show in figure [5] the zero-temperature phase diagram of the 2D spin-1/2 Ji- 
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Figure 4. The extrapolated CCM LSUBn results for the staggered magnetization, 
M, for J[ = 0.2,0.4,0.6,0.8,1.0 of the s = 1/2 Ji-J[-J2 model. The LSUBn results 
are extrapolated in the limit n ^ oo using the set n — {4, 6, 8, 10}. The NN exchange 
coupling Ji — 1. 




Figure 5. The extrapolated CCM LSUBn results for the gs phase diagram of the 
s = 1/2 Ji-J[-J2 model. The LSUBn results are extrapolated in the limit n — > oo 
using the set n = {4,6,8, 10}. The NN exchange coupling Ji = 1. QTP = quantum 
triple point. 

J[-J2 model, as obtained from our extrapolated results for both the gs energy and the 
order parameter. The phase diagram using the extrapolated LSUBn results based on 
n = {2, 4, 6, 8, 10} is very similar to that of figure which again adds credence to the 
validity and stability of our results. 

Within the high-order CCM that we have used our results certainly seem to provide 
clear and consistent evidence for a QTP at J[ ~ 0.6 for (Ji = 1). For J[ < 0.6 there 
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exist only the Neel and stripe phases, with a second-order transition between them. For 
J[ > 0.6 there exists an intermediate (magnetically disordered, i.e. paramagnetic) 
quantum phase, which requires further investigation. Although the nature of the 
intermediate phase is still under discussion, a valence-bond solid phase seems to be 
the most favoured [4, 9] . On the other hand, another possibility for the paramagnetic 
phase is the resonating valence bond (RVB) phase [16]. 

7. Discussion and conclusions 

In conclusion, our most important result is the evidence presented for a QTP for which 
our best estimate is J[ ~ 0.60 ± 0.03, J2 ~ 0.33 ± 0.02 for Ji = 1. Below this point we 
predict a second-order phase transition between the quantum Neel and stripe phases, 
whereas above it these two phases are separated by an intermediate phase. Although 
the order parameters for the Neel and stripe phases vanish on the Neel-intermediate and 
stripe-intermediate phase boundary lines respectively above the QTP, we are unable to 
conclude more about the nature of the transitions at those boundaries, since the present 
calculations have not addressed at all the nature of the intermediate phase. 

Other calculations on this model [16, 17] differ predominantly by giving a QTP at 
J[ = = J2 for Ji = 1. We believe that the difference arises essentially from the nature 
of the alternative methods used. For example, due to the small size of the lattices used, 
the ED calculations of Ref. [16] might miss the longer-range correlations that become 
increasingly important the nearer one approaches the QTP. By contrast, at any level 
of truncation the CCM always incorporates some long-range correlations through the 
important exponentiated parametrization of the wave function that lies at the heart of 
the method, as described in Sec. [3l Also, the results of Starykh and Balents [17] are 
based on analysis of a continuum version of the Ji-J[-J2 model rather than on the 
discrete lattice model itself. This could easily account for our differing predictions for 
the position of the QTP. 

As we have discussed above (and see figures H] and [5]), our results indicate that 
there is a continuous (second-order) transition below the QTP between the two ordered 
semiclassical phases with different broken symmetries. This is undoubtedly a rather 
novel and surprising result that has few precedents, and one that seemingly violates 
the Landau criterion of symmetry change. We note, however, that a similar scenario 
has been discussed recently in the related, but different, context of deconfined quantum 
criticality, as we explain more fully below. 

Thus, it has been argued recently by other authors [29] that there also exists a 
continuous (i.e., second-order) phase transition for the spin-1/2 pure J1-J2 model on the 
2D square lattice (i.e., our case with J[ = Ji) between the Neel state and what we have 
called here the intermediate state, and which those authors identify as a paramagnetic 
valence-bond solid (VBS) state. Such direct second-order quantum phase transitions 
between two states with different broken symmetries, and hence characterized by two 
seemingly independent order parameters is difficult to understand within the standard 
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Ginzburg-Landau critical theory, as we now discuss. 

Thus, the competition between two such distinct kinds of quantum order associated 
with different broken symmetries would lead generically in the Ginzburg-Landau 
scenario to only one of three possibilities : (i) a first-order transition between the two 
states, (ii) an intermediate region of co-existence between both phases with both kinds 
of order present, or (iii) a region of intermediate phase with neither of the orders of 
these two phases present. A direct second-order transition between states of different 
broken symmetries is only allowable within the standard Ginzburg-Landau critical 
theory if it arises by an accidental fine-tuning of the disparate order parameters to 
a multicritical point. Thus, for the spin-1/2 pure J1-J2 model and its quantum phase 
transition suggested by Senthil et al. [29], it would require the completely accidental 
coincidence (or near coincidence) of the point where the magnetic order parameter (i.e., 
the staggered magnetization) vanishes for the Neel phase with the point where the 
dimer order parameter vanishes for the VBS phase. Since each of these phases has 
a different broken symmetry (viz., spin-rotation symmetry for the Neel phase and the 
lattice symmetry for the VBS phase) , one would naively expect both that each transition 
is described by its own independent order parameter (i.e., the staggered magnetization 
for the Neel phase and the dimer order parameter for the VBS phase) and that the two 
transitions should therefore be independent of each other. 

By contrast, the "deconfined" type of quantum phase transition postulated by 
Senthil et al. [29] permits direct second-order quantum phase transitions between such 
states with different forms of broken symmetry. In their scenario the quantum critical 
points still separate phases characterized by order parameters of the conventional (i.e., 
in their language, "confining") kind, but their proposed new critical theory involves 
fractional degrees of freedom (viz., spinous for the spin-1/2 J1-J2 model on the 2D 
square lattice) that interact via an emergent gauge field. For the specific example of the 
spin-1/2 J1-J2 model the order parameters of both the Neel and VBS phases discussed 
above are represented in terms of the spinous, which themselves become "deconfined" 
exactly at the critical point. That the spinous are the fundamental constituents of both 
order parameters then offers a natural explanation for the direct second-order phase 
transition between two states of the system that otherwise seem very different on the 
basis of their broken symmetries. 

Despite the compelling nature of the arguments posited by Senthil et al. [29], we 
should mention, however, that other authors believe the phase transition in the J1-J2 
model not to be due to a deconfinement of spinous. For example, Sirker et al. [9] have 
given arguments, based both on numerical results from series expansion analyses and 
on spin- wave theory, that the spin-1/2 J1-J2 transition is not of the above second-order 
deconfined type but, in their view, is more likely to be a (weakly) first-order transition 
between the Neel phase and a VBS phase with columnar dimerization. 

One should also note that other, less radical, mechanisms have also been proposed 
to explain such second-order phase transitions and their seeming disagreement (except 
by accidental fine tuning) with Ginzburg-Landau theory. What seems clearly to be 
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minimally required is that the order parameters of the two phases with different broken 
symmetry should be related in some way. Thus, a Ginzburg-Landau-type theory can 
only be preserved if it contains additional terms in the effective theory that represent 
interactions between the two order parameters. For example, just such an effective 
theory has been proposed for the 2D spin-1/2 J1-J2 model on the square lattice by 
Sushkov et al. [30], and further discussed in Ref. [9]. 

Whether or not the deconfined phase transition theory of Senthil et al. [29] survives 
the controversy that still surrounds it, we note that one of the motivations that led to it 
was the existence of various other numerical calculations in recent years that also point 
to a direct second-order quantum phase transition between phases of different broken 
symmetry that are characterized by seemingly independent order parameters. Examples 
include the quantum phase transitions between : (i) an antiferromagnetic Mott insulator 
and a dx^-y^^. superconductor in a 2D Hubbard model on a square lattice, based on 
quantum Monte Carlo (QMC) simulations on lattices up to size 16 x 16 [31]; and (ii) 
superfluid and stripe-order phases in a 2D square-lattice spin-1/2 XY model with both 
a nearest-neighbour coupling term (J) and a four-spin ring exchange term {K), based 
on a QMC method (viz., a stochastic series expansion technique) with lattices up to size 
64 X 64 [32] . While most previous such numerical evidence for continuous (second-order) 
quantum phase transitions between states with different broken symmetry has come 
from QMC simulations, the well-known "sign problem" inherent to QMC techniques 
has meant that it has not been easy to apply that method to frustrated spin-lattice 
systems of the type considered here. We believe that the use of the CCM for such 
systems, as reported here, opens up a new arena and sheds a new spotlight on this 
fascinating and still unresolved larger field. 

Returning to our own results for the spin-1/2 Ji-J[-J2 model on the 2D square 
lattice, of course, one may also argue that what we have observed as a continuous 
(second-order) transition below the QTP may in reality be a very weak first-order 
transition, which is thereby not in violation of the Landau symmetry change criterion. 
Our completely independent calculations for the two semiclassical phases can, obviously, 
never entirely preclude this latter possibility. Nevertheless, it is clear from our results 
(and see figures [2] and Hj) that the data below the QTP are fully consistent only 
with a transition which if it is not continuous is very weakly first-order for all values 
< J( < 0.60 below the QTP. 

Finally, one could also argue that in this same range there might exist a very 
narrow strip of intermediate quantum-disordered phase, which would then reconcile 
our results (at least qualitatively) with those from the exact diagonalization of small 
clusters by Sindzingre [16]. While this possibihty, again, can never be entirely ruled out 
by any numerical calculation such as ours, we have clearly demonstrated that our own 
extrapolation schemes are both robust and internally consistent enough to rule out any 
but an extremely narrow strip of intermediate disordered phase for < J( < 0.60. 

We end by noting that two of the unique strengths of the CCM are its ability to deal 
with highly frustrated systems as easily as unfrustrated ones, and its use from the outset 
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of infinite lattices, which leads in turn to its ability to yield accurate phase boundaries 
even near a possible QTP. Our own results for the ground-state energy and staggered 
magnetization provide a set of independent checks that lead us to believe that we have 
a self-consistent and coherent description of this extremely challenging system. We 
also believe that the present high-order COM results are among the numerically most 
accurate for this and related spin-lattice models containing frustration. Nevertheless, our 
suggested novel phase scenario certainly needs further confirmation by the application 
of alternative high-order methods. It would also be of considerable interest to repeat 
the investigation for the computationally more challenging case of the same system for 
spin-1 particles, and we intend to report on this system in the future. 
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